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ABSTRACT 



> 

\q [ We calculate perturbatively the normalized spatial skewness, S3, and full 

three-point correlation function (3PCF), induced by gravitational instability 
O °f Gaussian primordial fluctuations for a biased tracer-mass distribution in 



flat and open cold-dark-matter (CDM) models. We take into account the 
explicit dependence on cosmological parameters, the shape and evolution 
of the CDM power spectrum, and allow the bias to be nonlinear and/or 
evolving in time, using an extension of Fry's (1996) bias-evolution model. We 
derive a scale-dependent, leading-order correction to the standard perturbative 
expression for S3 (Fry & Gaztanaga 1993) in the case of nonlinear biasing, 
as defined for the unsmoothed galaxy and dark-matter fields, and find that 
this correction becomes large when probing positive effective power-spectrum 
indices, i.e., scales above 100 h^ 1 Mpc for reasonable CDM models. This term 
implies that the inferred nonlinear-bias parameter, as usually defined in terms of 
the smoothed density fields, might in general depend on the chosen smoothing 
scale, and could allow better constraints on both the linear- and nonlinear-bias 
parameters on the basis of skewness measurements alone (or at least distinguish 
between the smoothed and unsmoothed bias pictures), if S 3 could be measured 
over very large scales. In general, we find that the dependence of S3 on the 
biasing scheme can substantially outweigh that on the adopted cosmology, 
with linear and nonlinear bias separately giving rise to distinct signatures in 
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the skewness, but degenerate ones in combination. We demonstrate that the 
normalized 3PCF, Q, is an ill-behaved quantity, and speculate that reported 
discrepancies between perturbative and TV-body predictions for Q may arise in 
part from systematic errors associated with the poor choice of normalization. 
To avoid this problem we investigate Qv, the variance-normalized 3PCF. The 
configuration dependence of Qv shows similarly strong sensitivities to the bias 
scheme as S3, but also exhibits significant dependence on the form of the CDM 
power spectrum. Though the degeneracy of S3 with respect to the cosmological 
parameters and constant linear- and nonlinear-bias parameters can be broken 
by the full configuration dependence of Qv, neither statistic can distinguish 
well between evolving and non-evolving bias scenarios, since an evolving bias is 
found to effectively mimic a smaller but constant bias. We show that this can 
be resolved, in principle, by considering the redshift dependence of ( which can 
also yield direct constraints on Qq and the epoch of galaxy formation. 

Subject headings: large scale structure of the universe — cosmology: theory - 
galaxies: clustering — galaxies: statistics 

1. INTRODUCTION 

Recovery of the primordial distribution of density perturbations is crucial to 
understanding the origin of large-scale structure. Inflation, possibly the most promising 
paradigm for the origin of structure, predicts Gaussian initial conditions (ICs), while 
most alternative models generally predict some deviation from Gaussianity. It is well 
known that gravitational instability will induce non- Gaussianity in an initially Gaussian 
distribution of cosmological perturbations, as the nonlinear collapse process gives rise to 
interactions between initially uncoupled modes and thus to non-zero higher-order moments. 
As fluctuations become nonlinear on progressively larger scales, larger objects become 
virialized, producing a hierarchy of structure. In the quasi-linear (QL) regime, where the 
rms fluctuations are small compared to unity, perturbative solutions to the fluid equations 
for a self-gravitating pressureless fluid in an expanding universe can be used to derive 
the resulting n-point correlation functions^ (CFs) and their Fourier-space counterparts, 
the polyspectra, induced by gravitational instability (Peebles 1980; Fry 1984; Goroff et 
al. 1986). In particular, the three-point correlation function (3PCF), or equivalently, its 
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Fourier transform (FT), the bispectrum, is the lowest-order intrinsically nonlinear statistic 
and can therefore place strong constraints on models of structure formation^. 

The perturbation-theory (PT) calculation predicts so-called hierarchical behavior of 
the polyspectra, whereby those of higher order can be simply related to the power spectrum, 
since the physics on QL scales is still determined entirely by the ICs. In the case of the 
bispectrum, 

B(kt, fc a , k 3 ) = Q(h, k 2 , k 3 ) [P(h)P(k 2 ) + P(k 2 )P(k 3 ) + P(k 3 )P(h)} , (1) 

where P(k) is the power spectrum and Q is of order unity and depends only on the 
magnitudes, ki —\ ki |, of the three wave vectors forming a triangle in /c-space. (Fry 1984; 
Goroff et al. 1986). A similar form was proposed, on purely empirical grounds, by Peebles 
(1974, 1975, 1980) for the spatial 3PCF, 

C(n2, r 23 , r 3 i) = Q(r 12 , r 23 , r 31 ) [£(ri 2 )6(r 23 ) + £(ri 2 )6(r 3 i) + £(r 23 )6(r 3 i)] , (2) 

where £(r) is the two-point correlation function (2PCF) and the normalized amplitude, 
Q, depends only on the physical separations, —\ — rj \ , of the three points in real 
space given by coordinates r i: and is assumed to be of order unity. It should be noted that 
though many authors assume Q to be a constant, this is strictly untrue, nor is Q given by 
the fc-space <2; rather it is a function defined by equation (|2]), independent of the validity of 
the hierarchical model. A related, more tractable statistic which is commonly used is the 
normalized spatial skewness, S 3 (R), related to Q evaluated at zero lag for a distribution 
smoothed over an effective scale half-width R. 

Studies of the spatial 3PCFs of different tracer populations on small scales (i.e., below 
10 h~ x Mpc) have generally yielded results consistent with the assumed hierarchical form 
above, i.e., values of Q ~ 1 ±0.5 (Groth & Peebles 1977; Jing & Zhang 1989; Toth, Hollosi, 
& Szalay 1989; Gott, Gao, & Park 1991; Jing, Mo, & Bdrner 1991; Baumgart & Fry 1991; 
Bouchet et al. 1993; Jing & Borner 1998) and values of 5* 3 of order a few (Gaztanaga 1992; 
Bouchet et al. 1993; Cappi & Maurogordato 1995; Gaztanaga, Croft, & Dalton 1995), and 
numerical simulations have also lent support to the hierarchical model (Bouchet, Schaeffer, 
& Davis 1991; Fry, Melott, & Shandarin 1993, 1995; Bernardeau 1994a; Baugh, Gaztanaga, 
& Efstathiou 1995; Colombi, Bouchet, & Hernquist 1996). At larger scales, however, the 
detailed behavior of these higher-order statistics, such as the configuration-dependence of 



4 Note that the n-point CFs are valid statistics only if the universe is homogeneous on large scales (as 
opposed to exhibiting, say, multi-fractal structure). We invoke the assumptions of large-scale homogeneity 
and isotropy, so that the CFs depend only relative distances and are independent of orientation. 
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Q, is not well constrained. Often, results for Q are quoted as constant simply because they 
reflect averages over all geometries, thus wiping out valuable information. Jing & Borner 
(1997, hereafter JB97) and Matsubara & Suto (1994) do investigate the shape dependence 
of Q, but find substantial disagreement between the predictions of QL PT and those of 
iV-body simulations. Jing & Borner (1998) further find significant devations between the 
observed configuration dependence of Q pro j, the projected Q perpendicular to the line of 
sight, and the standard hierarchical result. These differences may be attributable to various 
mechanisms, such as nonlinear effects, finite volume effects, and projection effects, which 
can reduce the configuration dependence of both observed and iV-body results for Q, as 
compared with QL PT predictions. The observed discrepancies, however, might also suggest 
that problems exist in the definition of Q. 

That the normalized amplitudes Q and S3 have been observed to exhibit roughly 
stable values near unity, as predicted, over scales on which ( and £ individually vary 
by many orders of magnitude, has been interpreted as loose support for the scenario of 
gravitational evolution of Gaussian ICs. It is clear, however, that the uncertainties in 
theoretical models, as well as in the data, are still appreciable. Much of this uncertainty is 
due to the fact that existing redshift surveys have not been both large and deep enough to 
yield strong constraints, particularly on QL scales. However, with the advent of new, deep 
surveys at various wavelengths (radio, optical, infra-red, x-ray), probing many different 
populations (galaxies, quasars, clusters, etc.), a rapidly growing database of spectroscopic 
and photometric redshifts, and powerful new computing resources, it will soon be possible 
to achieve precise measurements of higher-order spatial CFs. Interpreting the data from 
these upcoming surveys will require more realistic modeling of a variety of effects which 
may enter into the theoretical calculations. For example, the 3PCF can depend on the 
detailed form of the power spectrum (Frieman & Gaztahaga 1994; JB97) and, particularly 
for the case of angular CFs (see Buchalter, Kamionkowski, & Jaffe 1999, hereafter BKJ), 
on the density of nonrelativistic matter, the cosmological constant, A, or the density of 
some other form of matter (Bouchet et al. 1992; Bernardeau 1994a; Bouchet et al. 1995; 
Catelan et al. 1995; Martel 1995; Scoccimarro et al. 1998; Kamionkowski & Buchalter 
1998). Furthermore, smoothing of the density field, as occurs when determining S 3 from 
the moments of counts-in-cells, must be taken into account (Bernardeau 1994a,b). Perhaps 
most importantly, however, biasing of the tracer population will have a significant effect on 
the observed 3PCF (Fry 1984; Fry & Gaztahaga 1993, hereafter FG93). 

Bias — the notion that observed structures do not exactly trace the underlying density 
field — is advocated by both theory and observations, and may in general be scale and/or 
time dependent. In deterministic models, the fractional density perturbation 5(r) of the 
unsmoothed tracer- mass distribution may be expanded in terms of the perturbation 5 m {r) 
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to the total mass distribution, 

5(r)=M m (r) + |[5 m (r)] 2 + ---, (3) 

where bi is the linear bias term, 6 2 the first nonlinear term, etc., and 5 m is itself written 
as 5 m = 5$ + 8$ + • • • , where 8$ 8^" x \ 8$ is the linear solution, and 8$ is the 
leading-order departure from the Gaussian ICs. Note that the bias terms defined above are 
scale- independent, since equation @ is a local relation between the unsmoothed density 
fields at every point (we comment further on this assumption below). Like gravitational 
evolution or primordial skewness, the existence of nonlinear bias will also give rise to a 
non-zero third moment in the present-day density distribution of the tracer mass (since the 
3PCF is intrinsically second order, and thus requires the b 2 term in the PT calculation). 
These effects cannot be reliably distinguished by present data. Furthermore, though 
Gaussian ICs are sufficient to produce hierarchical behavior, they are not necessary; even 
if the data are shown to follow the predicted hierarchy to high accuracy, it is possible 
that this result is merely a coincidence of the bias. The impact of bias and its variation 
(thought to be primarily with redshift) on the 3PCF must therefore be taken into account 
if the predictions of Gaussian ICs are to be tested unambiguously. It has been shown that 
the detailed behavior of the 3PCF can constrain the initial density profile (Fry & Scherrer 
1994; Heavens 1998), distinguish between contributions to clustering from gravitational 
collapse vs. bias (Szalay 1988; FG93; Fry 1994; Frieman & Gaztanaga 1994) and, for a 
given tracer population, yield determinations of the linear and nonlinear bias parameters 
(Frieman & Gaztanaga 1994; Gaztanaga & Frieman 1994; Jing 1997; Matarrese, Verde, 
& Heavens 1997). Moreover, if different populations, such as clusters and galaxies, are 
differently biased relative to the underlying mass distribution, then measuring the 3PCFs 
of these populations can provide multiple, independent constraints. Several authors have 
recently argued that bias, and its evolution, is the dominant consideration in evaluating 
structure-formation models against clustering data, and have pointed to the need for 
employing more sophisticated bias models than the simplest ones often used (i.e., local, 
linear, and non-evolving) (JB97, Jing & Borner 1998; Bernardeau 1995; Matarrese et 
al. 1997). 

In this paper, we calculate the normalized spatial skewness, Ss(R), and normalized 
3PCF, Q(t*i, r-ii 1*3), in the QL regime assuming Gaussian ICs. We consider open and flat 
cosmological models with arbitrary allowed values of Qq and Q\, and take into account the 
dependence on the detailed form and evolution of the CDM power spectrum. We further 
include the effect of linear and nonlinear bias, and in particular, their time evolution, via 
an extension of the Fry (1996) bias-evolution model to the case of an arbitrary expansion 
history (BKJ perform similar calculations to obtain corresponding results for the angular 



- 6- 



skewness and 3PCF). For comparison, we derive the results for S 3 with a power-law 
spectrum, including a scale-dependent, leading-order correction to the standard expression 
for S3 (FG93) in the case of a nonlinear bias as defined in equation (|3|). This correction 
term becomes appreciable for positive effective spectral indices, corresponding to scales 
R > 100 h~ l Mpc for CDM models. The existence of this term implies that the value of the 
nonlinear-bias parameter, as defined for smoothed density fields, could in general depend 
on the adopted smoothing scale. For a given CDM model, this dependence, in principle, 
would allow a more accurate determination of the linear- and nonlinear-bias parameters on 
the basis of skewness measurements alone, provided however, that S 3 could be measured 
over sufficiently large scales. We show that Q is a poorly-defined statistic, exhibiting rapid 
variation and divergences which do not arise from the behavior of the 3PCF itself, but 
are rather due to the quantity [£(ri 2 )£(r 23 ) + £( r i2)6( r 3i) + £(r 23 )£(r 3 i)] in equation (g) 
acquiring values of or near zero in various cases [the fc-space Q defined in equation ([[[) 
is not subject to this normalization problem]. This dramatic behavior of Q is in marked 
contrast to that predicted by some iV-body simulations (JB97; Matsubara & Suto 1994) 
and we propose that this discrepancy is due at least in part to the practical inability 
to measure CFs to high precision near zero. We suggest that this may be addressed by 
considering iV-body results for Q at earlier epochs, where nonlinear effects can be safely 
ignored; to avoid this problem altogether, we define the variance-normalized 3PCF, Qy, 
given by ( divided by the square of the variance at a given scale. We find that S 3 and 
Qv do indeed depend strongly on the bias scheme, as has been suggested. In particular, 
a significant linear (nonlinear) bias term produces a relative decrease (increase) in both 
normalized amplitudes, and a significant linear term is found to reduce the variation of S3 
with smoothing scale and of Qv with triangle geometry, as compared with the unbiased 
predictions. While S3 is found to be only mildly sensitive to the cosmological parameters, 
Qv can vary appreciably with the value of the parameter r ~ QqIi, determined by the epoch 
of matter-radiation equality, that locates the peak in the power spectrum. In general, Qv 
is found to depend strongly on the triangle geometry, with smaller, elongated structures, 
corresponding to filaments or sheets, being more strongly clustered. The rich configuration 
dependence of the 3PCF is shown to break the degeneracy present in S 3 on scales R < 100 
h^ 1 Mpc, between the cosmological parameters and in particular the constant linear- and 
nonlinear-bias parameters. Both statistics, however, are unable to distinguish well between 
models of evolving and non-evolving bias, since an evolving bias effectively acts as a smaller, 
constant bias. We show that measurements of bias evolution, as well as of Qq and the epoch 
of galaxy formation, could potentially be isolated by considering ( as a function of redshift 
in a very deep (z > 1) survey. Alternatively, BKJ show that similar results can be obtained 
more feasibly from measurements of the angular 3PCF. 
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2. 



SPATIAL SKEWNESS 



The results below are based on leading-order PT and are thus restricted to the QL 
regime, where £ <^ 1. Though the leading-order contributions to higher moments of the 
density distribution are appreciable at these scales, higher-order nonlinear contributions 
are not expected to be significant, and various authors have shown that observations 
and simulations in this regime are well-matched by leading-order results alone (Szalay 
1988; Toth, Hollosi, & Szalay 1989; Gott, Gao, & Park 1991; Bouchet, Schaeffer, & 
Davis 1991; Fry, Melott, & Shandarin 1993, 1995; Colombi, Bouchet, & Hernquist 1996; 
Scoccimarro et al. 1998). In addition, numerical simulations suggest that large scales still 
obey leading-order PT even when small scales have become fully nonlinear (Bouchet & 
Hernquist 1992; Fry, Melott, & Shandarin 1995), and further that the results of QL PT hold 
even on scales where the rms fluctuation is ^ 1 (Bernardeau 1994a; Baugh, Gaztanaga, & 
Efstathiou 1995; Fry, Melott, & Shandarin 1995). 

We wish to derive an expression for the normalized 3PCF, Q, defined in equation (|2|), 
for a distribution of tracer masses, such as galaxies, clusters, or radio sources. This 
quantity has the advantage that it is independent of the overall normalization of the power 
spectrum. The full spatial 3PCF is given by C( r i ; r 2, fs) = (fi( r i)8( r 2) l 5( r 3)) , where r is the 
comoving position and 5(r) = [p(r,t ) — p(t )]/p(t ) is the present fractional perturbation 
to the unsmoothed tracer-mass density. Although the spatial 3PCF is nominally the 
function of nine quantities, statistical homogeneity and isotropy guarantee that the 3PCF 
depends only on three parameters which may be taken to be, e.g., the three distances 
between the three points. Still, evaluating Q over all allowed configurations may be fairly 
complicated, so it has become increasingly popular to focus instead on the normalized 
skewness S3 = (5 3 (r)) (5 2 (r)) 2 . This statistic is likewise independent of the power-spectrum 
normalization and is far easier to calculate than its full n-point counterpart, while still 
preserving much of the information contained therein. 

In practice, 6(r) cannot be evaluated continuously, since the observed tracer-mass 
distribution is discrete, so we define Sr(v) to be the density contrast smoothed over a cell 
with an effective comoving scale half-width In particular, we shall choose a spherical 
top-hat window function; in terms of the Fourier components, 5(k), of the unsmoothed 
density field, 




(4) 



5 We implicitly assume 5 = = 



1 oc £ — 1, where n is the number of discrete counts. 
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where W(x) = 3^7r/2 x^ 2 J 3/2 (x) is the FT of the three-dimensional spherical top-hat 
window function and J n (x) is a Bessel function of order n. Since we are working in the 
regime where (5 2 ) <C 1, fluctuations at or below the smoothing scale will not significantly 
affect our results (Bernardeau 1994a,b; see §3). 

The counts-in-cells spatial 2PCF is then given by 

Urn) = (fe(rOfe(r 2 )> = // ^ ^ e^ ^ k ^W(k lR )W(k 2 R) {^(k,)) . 

(5) 

The power spectrum of the tracer mass, P(k,t), is defined via 



[SikJSfa)) = (27r) 3 5 D (k 1 + k 2 )P(k, t), (6) 
where Sd is the Dirac delta function and k =\ k\ |. Putting this into equation (|5[), we have 

(SnirMr,)) = J ^ e« kr ^P(k, t)W 2 (kR), (7) 

where one integral vanishes under the requirement k\ = —k 2 , T\ 2 — V\ — r 2 , and we have 
made use of the the even symmetry of the window function. Taking r\ = r 2 = r, we obtain 
an expression for the variance, 

MO) = (Sl(r)) = J 0^P(k lt )W 2 (kR). (8) 

The linear solution to the equations of motion for the underlying density contrast has the 
separable form, 

€Hr,t) = D(t)5^(r,0), (9) 
where D(t) is the linear-theory growth factor, given (as a function of redshift) by 
5Q E(z) roo l + z < 



D ^ = ^r ± j z ^'[^jf E(z) = ] ji:m + ^ 1+Wi \ (10) 

where ^ Qi — 1 and Qi represents the contribution to the overall energy density from 
species % having equation of state p = Wipf\ for an Einstein de-Sitter universe, D(t) is 
simply the scale factor, a(t). One can then see from equation (Bf) that the spatial and time 



6 In this formulation, the curvature of the universe contributes an amount fix to the total energy density, 
and yields the term 0^(1 + z) 2 in the sum in equation (|To|); Thus, in a universe with, for example, a ratio of 
nonrelativistic-matter density to closure density of Slo and a cosmological-constant contribution to the total 
density of Cl A , E(z) = y/Sh(l + zf + fl A + (1 - ft - ^a)(1 + zf. 
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dependence of the power spectrum can likewise be separated. Assuming the unsmoothed 
tracer- mass density contrast to be related to the total-mass distribution via equation (^), 
we write the leading-order result for the power spectrum of the unsmoothed tracer-mass 
fluctuations as 

P(k,t) = AblD 2 (t)k n T 2 (k), (11) 

where A is the overall amplitude and T(k) is a model-dependent transfer function. 
Substituting equation ([11]) into (§), setting x = kR, and performing the angular 
integrations, we have 

UO) = (<£(«■)) = j~dxx^T\x/R)W\x). (12) 

We now follow a similar procedure for the counts-in-cells 3PCF, 

x (5(fc 1 )5(fc 2 )5(fc 3 ))iy(fci J R)iy(fc 2j R)iy(fc3 J R). (13) 



The bispectrum, B, of the unsmoothed tracer-mass distribution is defined via 

;*(fci)5(fc 2 )*(fc 3 )) = (27r) 3 fe(fc 1 + k 2 + fe 3 )S(A:i, fc 2 , fc 3 , *), (14) 



and is given, to leading order in PT, by (Fry 1984; Goroff et al. 1986; Matarrese, Verde, & 
Heavens 1997) 



B{h,k 2 ,k 3 ,t) = P(k u t)P(k 2 ,t)- 



1 + ^ + C o S e ( ^ + ^] + (1 - y) cos 2 9 
\k 2 h 



+|} + (cyc), (is) 

where 9 is the angle between ki and k 2 , and /x is a function of the expansion history 
(Bouchet et al. 1992; Bernardeau 1994a; Bouchet et al. 1995; Catelan et al. 1995; 
Martel 1995; Scoccimarro et al. 1998; Kamionkowski & Buchalter 1998). For an 
Einstein-de Sitter universe, \i = 3/7 (Peebles 1980). Bouchet et al. (1992) derive the 
approximation /x ~ (3/7)f2g 2//63 for an open universe. Kamionkowski & Buchalter (1998) 
find /x ~ (3/7)fi 1 ^ 14 ° f° r a universe with a cosmological constant [see also Bouchet 
et al. (1995)]. Bernardeau (1994a), Catelan et al. (1995), and Scoccimarro et al. (1998) 
give analytical arguments that further show that /x will depend only very weakly on the 
expansion history, and Kamionkowski & Buchalter (1998) verify this numerically for a 
variety of models. For plausible models with Qq > 0.1, the corrections to /x are below 2% 
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but we include them here for completeness. Note that the expression for B includes the 
dependence on the nonlinear bias term b 2 , since the first non-trivial terms in the connected 
n-point CFs require PT to order n — 1 [cf. equation (Q)]. Setting T\ = r 2 = r 3 = r, 
Xi = kiR, taking k\ to lie in the 9 = direction, and inserting equations flT4] ) and ( |ToD into 
(Pf), we obtain the skewness, 



QA b\D {t) r°o 



(2tt) 4 j R 2 ™+ 6 

x 



d6sm6{ — 
[bi 



roc roc 

/ dx l x1 +2 T 2 {x l /R)W{x x ) / dx 2 x 2 a+2 T 2 (x 2 /R)W{x 2 ) (16) 
Jo Jo 

(1 — /i) cos 2 9 



;i + /i) + cos#(^ + ^ 

X 2 Xi 



where a factor of 3 arises from symmetry considerations applied to the 2 permutations in 
equation (|15|). Noting that x 3 = \J x 2 + x\ + 1x\x 2 cos we can evaluate the 9 integral by 
using summation theorems for Bessel functions (Gradshteyn & Ryzhik 1980; Bernardeau 
1994a) which lead to the following relations: 



dO sin 3 9 W(x 3 ) 



d9 sin9 (1 + — cos#) W(x 3 ) 

X\ 



-W{x 1 )W{x 2 ) 
^X(x 2 )W( Xl ) 



(17) 
(18) 



- £(2j + 3/2)^ +3/2 (x 1 )^ + 3 /2 (x 2 ), (19) 

6 j=0 



where X(x) 
we have 



d9 sm9W(x 3 ) 

3sinx/x and W n (x) = 3yV/2 x~ 3/2 J n (x), so that W(x) = W 3 / 2 (x). Thus, 
l2A 2 b\D A (t) 



x 



(2vr) 4 J R 2 "+ 6 
( 1 



roo poo 

/ dx 1 x 7 l +2 T 2 {x l /R)W{x l ) / dx 2 x 2 l+2 T 2 (x 2 /R)W{x 2 ) 
Jo Jo 

2(1-//). 



^X(x 1 )W(x 2 ) + ^X(x 2 )W(x 1 



-W{x x )W{x 2 ) 



2 bo " 



+ o 7§ £(2j + 3/2)^ + 3/ 2 (^i)W 2j - +3 / 2 (^) . 

6 °1 j=0 / 



(20) 



For a given choice of b\ and 6 2 , the skewness obtained from equations ([121) and (|20f) 
effectively depends only on the smoothing radius, i? and the form of the power spectrum; 
the explicit dependence on fl and fl\ (through the dependence of /i) is very weak. 
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Fig. 1. — The function e(n) from equation (p3|). For —3 < n < the contribution of this 
term to S 3 in equation (^) is negligible, but increases rapidly for n > 0, due to the increased 
weighting of small-scale fluctuations by the unbounded power-law spectrum. 
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2.1. Effects of The Power Spectrum 

In general, we will consider a CDM transfer function given by (Bardeen et al. 1986) 

T() = ln(l + 2.34g)/(2.34g) 

[1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] 1/4 ' 

where q = k p /T, V ~ QqIi, and k p is the physical wavenumber in units of h Mpc -1 . The 
result for S3 is then obtained by using equation ([H]) together with (p0|) and (|l|). It is 
instructive to compare the CDM results with those obtained assuming a power-law power 
spectrum, P(k) oc k n . In this case, the x-integrals associated with the first two terms in 
parentheses in equation (EDI) can be evaluated analytically, and we obtain: 



S, 



(S 3 R (r)) . . 1 



34 6 f 7 



+ 3 | [!+£(„)], (22) 



where 



2 



00 



J£°dxX n 1 J3/2(x)J 2 j+3/2(x) 



e(n) = J E(2J + 3/2)^ ^ - — ' J . (23) 



3 i=i 



J °° dx X n - 1 J 3 /2(x)J 3 /2(x) 

Equations (22) and (^) are valid for n in the range —3 < n < 1. The individual terms 
in the sum in equation (Ej|) can be evaluated explicitly (Watson 1966), and we find that 
e(n) 1 for —3 < n < 0, but diverges for n > 0, as demonstrated in Figure [I]. Like the 
(3 + n) term, the e(n) term [which arises from the term proportional to 62 i n equation Q2DD] 
reflects the effect of the smoothing filter on the biased, contiguous field, defined in 

equation @. The filter is designed to separate smooth, large-scale fluctuations from the 
nonlinear, small-scale field. However, a nonlinear distribution having a scale-free power 
spectrum with n > contains so much small-scale power that the deviations induced by 
smoothing over these small-scale fluctuations become comparable to, or greater than, the 
perturbations above the smoothing scale, thus leading to the divergence in e(n). 



The form of equation (P2|) agrees with the previously known scale-free result for 5*3 
for a constant, nonlinear bias (Bouchet et al. 1992; Juszkiewicz, Bouchet, & Colombi 
1993; FG93; Bernardeau 1994a; Fry 1994), with the exception of the e(n) term. This 
difference arises from the fact that other authors calculating higher-order moments of the 
smoothed tracer- mass field (e.g., FG93) typically define the bias parameters using the form 
of equation ([|) to relate the smoothed tracer-mass and total- mass fields, i.e., 

S R = b i { ) 5 mR + ^-(5 mR ) 2 + ---, (24) 

where the b[ s ^ are the "smoothed bias" terms and a subscripted R again denotes a smoothed 
field. Biasing and smoothing do not commute; whereas the unsmoothed bias parameters 
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we define in equation (|3|) are truly local and scale- independent, the definition of bias in 
equation is inherently non-local, only fixing the smoothed bias parameters at the 
chosen smoothing scale. While it can be shown, using perturbation theory, that b± = b±, 
the higher-order smoothed bias terms will in general be different from their constant- valued 
counterparts in equation (§), and the magnitude of the difference will depend on the scale 
on which the b[ are defined. Had we used the smoothed bias parameters as defined in 
equation fl24|) , we would recover the result of FG93, i.e., equation (fffi) with b% — > b± , 
6 2 — > tyf , and e(n) = 0. Comparing our result for S3 with that of FG93, we can infer 
that &2 = fo 2 [l + £ ( n )}- In the no-smoothing limit, corresponding to n — > —3, we recover 
b^ — as expected. 

Of course, the unsmoothed bias parameters are defined only in the fictitious limit of 
continuous density fields. In practice, one might expect the efficiency of galaxy formation 
to depend on the matter density not only at the point of interest, but in the neighboring 
vicinity as well. Models such as peak biasing and halo biasing (Mo & White 1996) in fact 
imply a biasing relation such as equation (p4[), where the smoothed bias parameters do in 
fact tend towards constant values in the large-scale limit. Such a scenario would thus be 
inconsistent with the definition of the unsmoothed bias parameters in equation ([3D, and 
would yield precisely the FG93 result for S3, rather than equation fl2"2"|). Fortunately, these 
two situations can be distinguished on the basis of the large-scale behavior of S3, as we 
will illustrate. Unless otherwise stated, it is hereafter understood that bias shall refer to 
the relationship between the unsmoothed tracer-mass and dark-matter fields, as given by 
equation (tj). 



Naturally, equation fl2"2|) is /^-independent, as expected for a scale- free power spectrum. 
While the actual power spectrum for any tracer-mass distribution will not be simply a 
scale- free power law with n > 0, it may be wondered whether, in the unsmoothed biasing 
picture, the divergence seen in Figure [I] will affect the results for S3 on scales where the 
slope of the true (CDM) power spectrum is positive. Bernardeau (1994a) shows that the 
scale-dependent CDM result, without bias, can be recovered by using an effective index, 

rflog&(0) 

n R = j: 5 3, (25) 

a log it 

in place of n in the scale- free result for S^.f\ If we assume, for the sake of argument, that 
this holds true when calculating S3 assuming our nonlinear biasing model, we can illustrate 



By contrast, using an effective index denned by nu = [d\ogP(k)]/dlogk, evaluated at wavenumber 
k = c/R (where c is constant of proportionality chosen as to yield the best-fitting results), yields fractional 
errors in S3 of at least 5% (20%) for the smallest (largest) scales shown in Figure 
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Fig. 2. — The upper panel shows the results for ur(R) for the three cosmological models 
described (the low-f2o result is valid for both open models and flat models with a cosmological 
constant). The lower panel shows the variation of the e term in equation (|22"D with smoothing 
scale in these CDM models. Note that this term becomes comparable to unity on scales of a 
few hundred h~ x Mpc, where ur(R) > 0. The Figure should be interpreted only qualitatively, 
since using e[nn(R)} in equation (p2|) appears to overestimate the actual contribution of the 
higher-order terms in the sum in equation (BO) for CDM models. 
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the impact of the e(n) term as a function of smoothing scale for various CDM models. 
The upper panel of Figure shows n R (R) for for three flat (Q + Q\ = 1) cosmological 
models: standard cold dark matter (SCDM; Q = 1, h = 0.5, n — 1), a tilted CDM model 
(TCDM; n = 1, h = 0.5, n = 0.8), and CDM with a cosmological constant (ACDM; 
Qo = 0.4, h = 0.5, n = 1), while the lower panel shows e[n R (R)] for these models. One 
can see from equations (|T2|) , (f23|), and (|25| ) that the results from the ACDM model will be 
identical to those for an open CDM model with the same value of Qq (OCDM; Q = 0.4, 
h = 0.5, n = 1). Though it appears that using ur in the e term in equation (|22|) , with 
our nonlinear-bias model, is not strictly correct for scale-dependent power spectra (see 
Figure |5|), we can at least qualitatively infer from Figure [| that the contribution of the 
e term to S3 can become significant on scales R > 100 h~ x Mpc, particularly for models 
with higher values of T. Thus, when considering data on the moments of counts in cells, 
smoothed on these scales, this correction term must be included to properly account for 
the presence of nonlinear bias. BKJ obtain a similar correction to the angular skewness, 
and corresponding corrective terms should arise in leading-order perturbative calculation of 
higher-order moments, such as the kurtosis, S4, and so on. For smoothed bias, as defined in 
equation (0), this result implies that if and b 2 are defined at some smoothing radius 
Ri, then at a sufficiently large scale R2 7^ Ri, b 2 will acquire a different value due to the e 
induced by the change in R. For a given CDM model, measurements of S3 using different 
smoothing radii should thus yield identical results for b\ = bi and b 2 , but different results 

(s) RI 

for b 2 .(] This suggests that the linear and nonlinear bias terms may be distinguishable, in a 
model-dependent fashion, on the basis of skewness measurements alone, provided however, 
that S3 can be measured over scales sufficiently larger than 100 h~ l Mpc. An illustration 
of this point is shown below, in Figure f| If, however, measurements do not show any 
disagreement between 6 2 and b^ on large scales, this would provide evidence to support 
scenarios such as peak-biasing or halo-biasing models and argue against biasing as defined 
in equation @. 

Since there is no known method of obtaining analytic results for the general case of the 
full 3PCF, assuming scale-dependent CDM power spectra, we will hereafter, for consistency, 
only present results obtained by numerical integration, unless otherwise noted. Many of 
the integrals throughout this paper involve highly oscillatory integrands, making them 
difficult to evaluate. We have checked that our numerical results for S3 in particular agree 
to high precision with those obtained using Bernardeau's (1994a) analytic prescription, 
and that our results in general are reliable to within a few percent. Figure |3] displays our 
numerical results for S^R) for an unbiased (pi — 1, b 2 — 0) tracer-mass population in four 



'We thank Roman Scoccimarro for clarification of this point. 
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Fig. 3. — The predicted normalized skewness, Ss{R), for an unbiased tracer-mass distribution 
smoothed with a spherical top-hat filter of radius R, for the four flat and one open 
cosmological models described in the text. The scale dependence of S3 arises from the 
adopted CDM transfer function. By comparison, the thin horizontal lines are the results 
obtained using power-law spectra with n = — 2 and n = —1.2. Note that the dependence of 
S3 on the assumed cosmological model is fairly weak, particularly if the observed constraint, 
0.2 < il h < 0.3, is imposed. 



-17- 



flat (J] +Q\=1) cosmological models: SCDM, TCDM, ACDM (each from above), and 
CDM with a high Hubble parameter (HCDM; Q = 1, h = 0.75, n = 1), as well as the open 
CDM model (OCDM, from above). 

As expected, Figure |3] shows that the CDM power spectra introduce a dependence on 



the smoothing scale, R, and [compared to equation (|22|) 1 a more substantial dependence 
on the combination of cosmological parameters in T. Still, the different cosmological 
models represented, though spanning a fairly broad range, yield fairly similar results for S3, 
certainly all consistent within the typical errors of current observations (e.g., Gaztanaga 
1992; Jing & Borner 1998). The slight variation seen is due primarily to the dependence on 
T; there is little difference between the SCDM and TCDM models (which differ only in the 
value of n) and virtually no difference between the OCDM and ACDM models, as expected 
from the very weak \i dependence. If the analysis is restricted to the range 0.2 < T < 0.3, 
as suggested by current data (e.g., Bartlett et al. 1998), one can infer from the Figure that 
the predicted values for S3 are relatively insensitive to the adopted cosmological model. 
In the R — ► limit, the untilted (n = 1) CDM curves tend toward the expected value of 
(34/7) = (34/7) — (3 + n k ) obtained with n k = \d\ogP{k)]/ dlogk = —3, the value of the 
effective spectral index of the (untilted) linear-theory CDM power spectrum at small scales; 
for tilted models, S3 tends toward a value n k — 1 smaller. At large R, they tend toward 
(6/7) = (34/7) — (3 + rife), where n k = 1. These perturbative results are expected to be valid 
at scales R ^ 5 — 10 h^ 1 Mpc, although as mentioned above, there are reasons to expect 
them to be reliable even at smaller scales. In practice, a comparison of these predictions 
with data will be limited at small scales by increasingly significant nonlinear effects, as 
well as Poisson shot noise arising from the discreteness of the counts-in-cells (Peebles 1980, 
Gaztanaga 1994), and at very large scales by sampling noise arising from small number of 
independent cells on the sky. Note for comparison, that the thin horizontal lines, which 
show the (virtually f2 -independent) analytic results for scale-free power spectra with 
n = —2 and n = —1.2 [the canonical index value obtained from measurements of the 2PCF 
at small scales (Peebles 1980)], provide poor fits to the CDM predictions. 



2.2. Effects of Bias and its Evolution 

The calculations above can be generalized by allowing the bias parameters to evolve 
with time, as advocated by theory and observations. Peacock (1997) predicts that the linear 
galaxy-mass bias term evolves from a value of roughly 6 at a galaxy-formation redshift of 
Zf ~ Q — 8, to a nearly unbiased value of ^ 1 today, with all models of structure formation 
requiring some degree of bias at z ^ 3. Matarrese et al. (1997) explore several different 
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models of bias evolution, and also find that the bias for objects of between 10 10 and 10 13 
M Q evolves from a value of a few at z ~ 5 to roughly unity today. Studies of the clustering 
strength of tracer populations that span an adequately large range of redshifts, such as radio 
galaxies and Lyman-break galaxies, show evidence for a large bias at high redshift, which 
decreases with time in a manner consistent with the predictions above. (Steidel et al. 1998; 
Cress & Kamionkowski 1998). Assuming that objects form at a fixed redshift Zf by some 
arbitrary local process which induces a bias at that epoch, and are subsequently governed 
purely by gravity, Fry (1996) derives the result for the bispectrum of a tracer population 
in an Einstein-de Sitter universe, including the time-dependence of the bias parameters. 
Matarrese et al. (1997), refer to this as the "object-conserving" model, since it does not 
account for merging (see below). Tegmark & Peebles (1998) generalize the Fry model to 
the case where the dimensionless mass-galaxy correlation coefficient may be different than 
unity, but we take the Fry model as a reasonable starting point for the exploration of the 
effects of bias evolution. 

It is straightforward to generalize the bias-evolution model of Fry (1996) to an arbitrary 
expansion history. Doing so, we obtain for the bispectrum in these models, 



Ci(*) + c 2 (t) cose + ^) + c 3 (t) cos 2 e 



B(k u k 2 ,k 3 ,t) = P(k u t)P(k 2 ,t) 

+2 permutations, (26) 



where 



r m (10/7)d 2 (t) + 2(6 1 .-l)[d(t)-2/7] + 6 2< 

Cl(t) ~ TO + ^-iP ' (27) 

C (t) (4/7)TO + 6 u -l] 

C3{t) ~ [d(t) + K-i]* ' (29) 

d(t) = D(t)/D(t*), and a subscripted asterisk denotes the value of that parameter at the 
epoch of formation. Note that we have now ignored the very weak dependence of the 
bispectrum on //; the dependence on the expansion history, i.e., on the species contributing 
to the total energy density, is contained in the growth factor, D{t). 

In this model, the linear-bias term effectively evolves as 

decaying towards unity, as in most bias models, as the observed matter settles into the 
potential wells of the underlying distribution. We note that if, in reality, merging of the 
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tracer-mass population does in fact play a considerable role, this could result in an anti-bias 
(61 < 1). Anti-bias has been seen in some TV-body simulations (Jenkins et al. 1997), and 
certain models of structure formation do require some degree of anti-bias to reconcile with 
current observations. The above model for bias evolution, however, will never yield b\(t) < 1 
if bi t > 1. To explore this possibility, one might wish to employ other bias-evolution models 
which account for merging. Here, however, we are interested primarily in the qualitative 
dependence of the results on the dominant form of bias variation, namely, its decay with 
time from initially large values at high redshift. 

The Einstein-de Sitter result for S3 with bias evolution, assuming a scale-free P(k) oc k n 
is presented in Fry (1996). For arbitrary E(z), the result becomes 

_ ^ 2 (t)(34/7) + (K - l)[6rf(t) - 8/7] + 3fe 2 , d(t) 
3 " W) + b u ~ I] 2 {3 + n) d(t) + b u -l (31) 

Here, we extend the calculation to the case of true, scale-dependent CDM power spectra, by 
using the bispectrum of equation fl26|) in (|17D, which amounts to making the substitutions 



>C7 2 (t), V^^^W, 3— C7 1 (t)-2C7 8 (t) + C7 8 (t) l (32) 

01 Oi o x 

in equation (p0|) . 



Figure [| displays the results for five different bias scenarios, each employing the SCDM, 
ACDM, and OCDM models above: an unbiased scenario (pi = 1, &2/&1 = 0, taken from 
Figure 2), a non-evolving, linear bias (b\ = 2, &2/&1 = 0), a non-evolving, nonlinear bias 
(61 = 2, b 2 /b\ = 1), an evolving, linear bias (&i„ = 5, &2*/&i = 0), and an evolving, nonlinear 
bias (fei, = 5, &2*/&i„ — 1)- For the non-evolving bias scenarios, the ACDM and OCDM 
models are taken to be identical, as suggested by Figure 2. These models do, however, differ 
in the case of evolving bias due to the different time dependence of the linear growth factor, 
and the differences between them would increase with decreasing Qq. For the evolving cases 
we assume a formation redshift of Zf — 5. In calculating the expressions involving 62, we 
must evaluate sums of the kind in equation (|20|). Numerical tests show that these converge 
rapidly, typically requiring at most three terms to achieve percent-level accuracy. 

These models are intended as an illustrative, but by no means exhaustive, representation 
of the possibilities. Even so, it is clear from Figure ^ that that S^(R) depends far more 
sensitively on the overall form of the bias than on the values of T or Qq within an individual 
biasing scheme. For example, adding a very slight linear bias to the low-r model would 
yield a result similar to the unbiased SCDM prediction, within current error estimates. 
Generally, the presence of any significant linear bias, b\ > 1, flattens the 6*3 (i?) curve, as 
compared with the unbiased case (solid lines). Furthermore, an observed Ss(R) curve which 
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Fig. 4. — The predicted normalized skewness for the five different bias scenarios described in 
the text. For each case, the SCDM, ACDM, and OCDM models are shown in thick, medium, 
and thin lines, respectively (for the non-evolving bias scenarios, the latter two models are 
identical). Note that the dependence on the overall biasing scheme far outweighs that on 
the cosmological parameters within a given scheme. In particular, a linear (nonlinear) bias 
term lowers (raises) the predicted curve with respect to the unbiased case, and a significant 
linear bias term flattens the curve appreciably. 
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is well below the predicted unbiased result can only be achieved by a significant linear bias 
term (assuming b 2 > 0), while one well above this value can only arise from the existence of 
a non-zero b 2 term or from anti-biasing (61 < 1). In general, the curves for an evolving bias 
produce less drastic shifts than their non-evolving counterparts, despite the fact that the 
respective bias terms are initially larger. This is because the evolution towards an unbiased 
state effectively mimics a smaller, constant bias [c.f. equation fl30l) 1. Note that the SCDM 
models with nonlinear bias exhibit an upturn at large R, arising from the contribution of 
the higher-order terms in the sum in equation (|20|) . The other CDM models with nonlinear 
bias show an upturn only at larger scales, as suggested by Figure 0. 

Figure [| also reveals that, despite the separate dependences of S3 on b%, b 2 , and their 
evolution, different biasing scenarios may be very difficult to distinguish on the basis of 
the skewness measurements alone. Different sets of parameter choices could easily yield 
curves for the various scenarios which essentially overlap, requiring high-accuracy data to 
differentiate them. An example of this is shown in Figure [5], where we see that, for scales 
R < 100 h^ 1 Mpc, the results from an SCDM model with constant, nonlinear bias, an 
OCDM model with evolving, linear bias, and an OCDM/ACDM model with constant, linear 
bias, are very similar. The OCDM model with evolving, linear bias is taken from Figure |], 
and according to equation (|30"D, is expected to be similar to an open model with a constant 
value of b\ = 2.04, or a ACDM model with a constant b\ = 1.8, and indeed Figure |5] shows 
this to be the case. For an inferior data set, this degeneracy between the cosmological 
parameters, and in particular b\, b 2 , and their time dependence, cannot be broken with 
the skewness, measured on scales R < 100 h" 1 Mpc, as it provides only one constraint on 
the combination of these quantities. Note, however, that for R > 100 h~ l Mpc, the SCDM 
model with nonlinear bias can be better distinguished due to the different scale dependence 
arising from the higher-order terms in the sum in equation (|20"D For comparison, the thin 
dot-dashed line shows the same result obtained using e[nji(R)} in equation (p2|), which 
appears to overestimate the numerical result (thick lines) obtained by evaluating the the 
sum in equation fl2"0"|). If this large-scale variation induced by nonlinear biasing is observed, 
it would allow better constraints on both the linear and nonlinear bias parameters solely 
from measurements of S3. The absence of this behavior would still provide valuable insight, 
arguing against our bias model and in favor of smoothed-bias models, where b\ and b 2 
tend towards constants at large scales. However, it is at present unlikely that 5*3 could be 
measured on such large scales with sufficient precision to address these issues. 
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Fig. 5. — An example of nearly degenerate predictions for S^R), arising from the fact that 
the skewness provides only one constraint on the combination of the cosmological and bias 
parameters, and the evolution of the latter. In this case, the SCDM model with constant 
nonlinear bias, OCDM model with evolving linear bias, and OCDM/ACDM model with 
constant linear bias are nearly identical out to smoothing scales of R ~ 100 h^ 1 Mpc. Above 
this scale, the SCDM model with nonlinear bias can be better distinguished due to the 



different scale dependence introduced by higher-order terms in the sum in equation ( 20f) . 
For comparison, the thin dot-dashed line shows the SCDM result using equation ( p2| ) with 
e[nfl(i?)], and appears to overestimate the impact of the nonlinear term at large scales. 
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3. SPATIAL THREE-POINT CORRELATION FUNCTION 

The one-point statistic S3, is a volume-averaged quantity and thus, while it is seen 
to preserve information about the overall dependence of clustering strength with scale, 
discards detailed information about the configuration dependence contained in the full 
3PCF. It has been shown, for example, that the configuration dependence of the bispectrum 
and 3PCF in leading-order PT can be used to separate the contributions of gravitational 
clustering and bias to the observed power (FG93; Fry 1994; Frieman & Gaztanaga 1994), 
as well as obtain multiple, independent constraints on the combination of the constant b\ 
and 62 parameters (Jing & Zhang 1989; Gaztanaga & Frieman 1994; Fry 1994, 1996; Jing 
1997; Matarrese, Verde, & Heavens 1997). 

Much attention has therefore been focused on the normalized spatial 3PCF, 

n( \ _ C( r 12, ?~23, r 3l) /qq\ 

g(ri2! r23 ' r3l) " arn)ar,s) + arunr,) + a^n^y (33) 

where the lengths form a triangle in real space. We can evaluate this expression using 



equations (||) and (jIBD , where we no longer include the window functions since we presume 
that the CFs are obtained from direct counting, rather than from smoothed counts in 
cells. In this case, it is more natural to use the unsmoothed bias parameters defined 
by equation ([|), rather than the smoothed bias parameters of equation (24). Using the 
techniques and conventions of the previous sections, we obtain, after some manipulations, 
the spatial 2PCF, 9 

/-I \ 1 f JJ j2r,n\ S ^ nkr Ab\D 2 (t) f n+2 m2f I \ sin X tn a\ 

i{r) = 2^J dkk P{k) — = 27rV*3 J dxx n + 2 T 2 (x/r) — , (34) 

and 3PCF, 

A 2 h^D 4 (t) 1 r°° r°° 
C(ri 2 , r 23 ,r 31 ) = \ 1 > n+3 n+3 / dx 1 x^ 2 T 2 (x l /r 31 ) dx 2 x n 2 + 2 T 2 (x 2 /r 23 ) 

4t7T ^*31 ^23 



9 Had we wished to calculate the 2PCF using counts in cells, there would be an additional factor 
of W 2 (xR/r) in the integrand, and corresponding terms in the 3PCF, which would alter the results 
on scales where R/r is non-negligible (i.e., separations comparable to the cell size), to account for the 
effects of smoothing on correlations at these scales. The dependence of the CFs on the ratio between 
the smoothing half-width and the pair separation can be appreciable; for scale-free power spectra with 
n = -2.2, -1.2, -0.2, 0.8 this term adds corrections to £ of 1%, 10%, 49%, and 294% for r = 2R, dropping 
to 0.14%, 1.2%, 4.3%, and 2.9% for r = 5R. For a CDM power spectrum, with T = 0.5, the corrections 
are more significant: 42%, 11%, and 2.4% for r/R values of 2, 5, and 10, respectively. For counts-in-cells 
CFs, one must take care to restrict measurements to scales significantly larger than the smoothing radius, 
or include the smoothing terms in the theoretical calculations. 
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Fig. 6. — The variation of triangle geometry with u and v, for a fixed value of r given by 
the side of the u — 1, v — equilateral triangle. Triangles with u — 1 are isosceles with the 
shorter sides equal, those with v = are isosceles with the longer sides equal, and v — 1 yields 
a colinear configuration. Larger values of u and/or v pick out more elongated structures. 
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1 + fi b 2 \ sinxi sinx 2 



b\ b\) x\ x 2 

cosfl 12 (—x\V{x 1 )V{x 2 ) + — x 2 2 V(x 1 )V(x 2 )) 
Vr 3 i r 23 * / 



cos 2 9 12 Y( Xl )Y(x 2 ) + 3V( Xl )V(x 2 ) 

+ (cyc), (35) 



6i 

-y(xi)y(x 2 ) - y(x 2 )y(x! 



where cos6 l i 2 is the angle between r i3 and r 23 , V(ar) = W(x)/3, and and 
= W(x) — (smx)/x. For an evolving bias, we simply change 

C 2 (t), If£— C s ( t ). (36) 




h 

Equation (|3^) does not contain terms analogous to the higher-order terms in the sum 
in equation (p0|) [or to the e{n) term in equation ([221)1, since evaluation of the 3PCF by 



direct counting requires no smoothing. With our bias model, such terms would be induced, 
however, from the calculation of the 3PCF using counts in cells (see footnote). One must 
be careful to realize that the nonlinear-bias parameter measured from the direct-counting 
3PCF in equation fl3~5| ) will be different from the smoothed value inferred from the FG93 
result for S3 (or, equivalently, from that inferred from the counts- in-cells 3PCF), though 



the two can be related, e.g., via the e(n) term in equation (25). 

The dependence of £(r) on r~( n+3 ) in the prefactor of equation (0) is the familiar 
result, yielding £(r) oc r~ L8 for a scale-free power spectrum with the canonical 
value of n = —1.2 (Peebles 1980). For more realistic CDM power spectra, there 
is an additional dependence on r arising from the scale dependence of the transfer 
function. Similarly, had we ignored the shape- dependent terms in the bispectrum 
[equation (pT5[)1, we would arrive at the standard empirical "approximation" for the 3PCF, 
C(r 12 ,r 23 ,r 31 ) oc £(ri 2 )£(r 23 ) + €(r 12 )€(r 31 ) + f (r 23 )£(r 3 i). However, the true scale and shape 
dependences of the 3PCF, as manifested in equation fl35D , can be exploited to gain far more 
information (Fry 1984; Jing & Borner 1997). Following Peebles & Groth (1975), we take 

^23 ^31 — r 23 , on \ 

r = r 12 , u = — , v = (37) 

n 2 r 12 

as the defining parameters for a triangle, where r\ 2 < r 23 < r 3 i, so that u > 1 and < v < 1. 
For this convenient choice of parameters, r fixes the overall size of the triangle, while u and 
v determine the exact size and shape. Figure [] shows the variation of the geometry with u 
and v for a given choice of r. 
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Fig. 7.— Predicted normalized 3PCF, Q(r,u,v), for the T = 0.5 (SCDM, solid lines) and 
T = 0.2 (OCDM/ACDM, dashed lines) models without bias. The columns represent values 
of r = 4, 8, 16, and 32 h^ 1 Mpc (left to right) while the rows represent values of u — 1, 
2, 3, and 4 (top to bottom). Note the different scalings on the vertical axes. Over these 
scales, the SCDM model shows dramatic variation of Q with triangle geometry, exhibiting 
divergences in some cases. 
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3.1. Results with no Biasing 

3.1.1. The Problem With Q 

Figure |7| shows the normalized spatial 3PCF, Q, plotted as a function of r, u, and 
v, using the transfer function of equation (^l|) for the T = 0.5 (SCDM) and T = 0.2 
(OCDM/ACDM) models, assuming no bias. The weak \x dependence makes the OCDM 
and ACDM models nearly indistinguishable on the scale of Figure [7], so only the former 
is plotted. Note the different scalings and the dramatic, non-monotonic variation and 
divergences seen for certain configurations, particularly in the SCDM model. JB97 also 
derive QL PT predictions for Q in CDM models and find very similar results for its 
configuration dependence, though they do not show cases where the behavior is particularly 
egregious. Comparing these results with those from iV-body simulations, however, JB97 
and Matsubara & Suto (1994) find that the predicted differences between various CDM 
models in iV-body simulations are markedly less than those predicted by QL PT, and 
moreover that in simulations, these models in general show far less rapid variation of Q 
with geometry than is seen in the QL PT predictions. JB97 claim that this discrepancy is 
due primarily to the absence of nonlinear corrections in the PT calculation, but we suggest 
another possibility. In the top panel of Figure § we show our prediction for Q(v) in a 
T = 0.5 model for r = 8 h~ x Mpc and u — 4. The lower panel reveals that the divergence in 
Q is caused not by any irregular behavior of ( (which is seen to vary in a smooth, monotonic 
fashion), but rather by the behavior of the denominator in equation (|3"3"|), which vanishes 
in this case at v ~ 0.325 (note that we have multiplied the numerator and denominator 
of equation (^) by 10 10 and 10 12 , respectively, for clarity). In this sense, Q is a poorly 
defined quantity, whose variation may not reflect any rapid variations in (, but rather 
stems from the manner in which ( is normalized. For the cases shown in Figure [7|, Q 
is particularly rapidly varying in the SCDM model only because the region of the power 
spectrum effectively sampled by a value of V — 0.5 happens to yield vanishing behavior 
of ^(^12)^(^23) +^(^12)^(^31) + £( r 23)£(?"3i) for some of the chosen configurations; similar 
behavior would be seen in other models at appropriately different scales. 

We speculate that the above-mentioned disagreement with iV-body predictions, which 
show far less variation in Q, may simply be due to the practical inability to measure 
£(^12) £(^23) + 6( r i2)£( r 3l) +^(^23)^(^31) to high precision near zero, thus smoothing out the 
iV-body result for Q. This conclusion is supported by the fact that all previous investigations 
find that the QL PT predictions for S3, which is normalized by a positive-definite quantity 
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Fig. 8. — The top panel shows the predicted Q(v) for r = 8 h^ 1 Mpc and u = 4, while the 
bottom panel shows the separate behaviors of the numerator (multiplied by 10 10 ; dashed 
line) and denominator (multiplied by 10 12 ; solid line) of equation ( pBf ) for this configuration. 
The 3PCF itself is smoothly varying; the large values of Q arise from the small values of the 
denominator, which passes through zero at v ~ 0.325. 
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(i.e., the square of the variance), do agree with iV-body results over the scales of interest^ 
(Bouchet, Schaeffer, & Davis 1991; Fry, Melott, & Shandarin 1993, 1995; Bernardeau 
1994a; Baugh, Gaztanaga, & Efstathiou 1995; Colombi, Bouchet, & Hernquist 1996; JB97). 
We propose that this explanation for the discrepancy in Q can be evaluated against the 
nonlinear-effects argument of JB97 by considering iV-body results at earlier epochs, where 
nonlinearities in the distribution of fluctuations are known to be small. We point out that 
the Fourier-space Q, defined in equation fll]), is normalized by a positive-definite quantity, 
and does not suffer from this problem entailed in the definition of the real-space Q. 



3.1.2. Results for the 3PCF 

To avoid the above difficulties altogether, we define Qy = (/[£,8h- 1 M P c(fy} 2 'i i- e - ; we 
normalize the 3PCF (like S3) by the square of the variance at a fixed smoothing scale, 



calculated for each model using equation (0). Thus, Qy is still independent of the 



power-spectrum normalization, but does not exhibit the rapidly-varying or divergent 
behavior seen with Q, since we are now dividing by a positive-definite quantity. The 
adopted smoothing scale is arbitrary; 8/i -1 Mpc is chosen merely for convenience. 

Figure g illustrates our predictions for Q v (r,u,v) for the SCDM, OCDM/ACDM, 
and TCDM models with no bias. Note that Qy is everywhere a smoothly-varying and 
well-behaved function. Since, in each model, Qy is normalized by a fixed quantity, the 
Figure clearly reveals the expected decrease in the 3PCF with increasing r. This behavior, 
however, is not simply the r _2(ra+3 ) variation predicted by the hierarchical model, which 
also fails to produce the shape dependence of the 3PCF. In general, we see that, within 
a given cosmological model, Qy increases with increasing v. Physically, this implies that 
clustering in the quasi-linear regime favors colinear structures (see Figure []), which can 
be identified, in the early stages of structure formation, with collapsed structures such as 
sheets and filaments. The decrease in Qy with increasing u is due to a combination of the 
scale and shape dependences; larger values of u pick out both larger and more elongated 
structures, giving rise to competing effects. Under our chosen normalization convention, 
the former effect wins out. Fry (1994, 1996) finds similar results by considering the fc-space 
Q, defined in equation ([]]). Though the CDM models shown all exhibit these common 
trends, the precise behaviors predicted for Qy by models with different values of T are 
quite distinct, owing to the different locations of the peaks in their power spectra. Overall, 



10 Hui & Gaztanaga (1998) note a separate measurement bias, applicable to S3 as well, arising from the 
fact that the ratio of two unbiased estimators is not unbiased estimator. 
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Fig. 9.— Predicted Q v (r,u,v) for the unbiased SCDM (solid lines), OCDM/ACDM (short- 
dashed lines), and TCDM (long-dashed lines) models. The columns represent values of r = 4, 
8, 16, and 32 hr x Mpc (left to right) while the rows represent values of u — 1, 2, 3, and 4 
(top to bottom). Note the different scalings on the vertical axes. All three models indicate 
that clustering strength decreases with scale and increases with increasing v. Models with 
different values of T, however, yield distinct predictions for Q v , though the SCDM and 
TCDM models, which differ only in the value of n, are not well distinguished on large scales. 
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Figure ^ suggests that measurements of the 3PCF on these scales, though not very sensitive 
to changes in the value of n (see BKJ), may discriminate well between high-r and low-r 
models. 



3.2. Effects of Bias and Its Evolution 

Figure |10| shows Qv{v) for the five bias scenarios from §2.2 for the OCDM and ACDM 
models (taken to be identical in the non-evolving bias cases) for r = 8 h^ 1 Mpc and 
u = 4. This represents only a small region of the parameter space we are exploring, but 
we find similar dependences of the 3PCF on the bias scheme as were seen for S3, namely a 
general flattening and reduction of Qv with increasing linear bias, and a relative increase 
with increasing nonlinear bias as compared with the corresponding linear scheme. In the 
evolving cases, the difference between these two models, which arises from their different 
linear growth factors, would be greater for lower values of Qq. 

Whereas the skewness gave only a single relationship, Ss(R), with which to distinguish 
different cosmological and biasing schemes, the full 3PCF allows much stronger constraints 
from consideration of the geometric dependence [i.e., a graph such as Figure |X0| for every 
pair of values (r, u)]. We illustrate this in Figure [11], where Qy{r,u,v) is plotted for the 
same three (nearly degenerate) models depicted in Figure |5]. The constant linear- and 
nonlinear-bias schemes (dotted and dot-dashed curves, respectively), which gave nearly 
identical results for 5*3, yield distinct predictions for Q v , in part because of their different 
values of T as well. It is well-known that measurements of the bispectrum and 3PCF can 
be used to obtain independent constraints on the constant-valued linear and nonlinear bias 
parameters (FG93; Fry 1994, 1996; Jing 1997; Matarrese, Verde, & Heavens 1997); this 
can be seen here from the different manners in which these terms appear in equation (|3"5|). 
The low-r models in Figure |TT|, with linear, evolving (dashed curves) and non-evolving 
(dashed curves) bias, however, remain fairly similar, since an initially large, evolving value 
of b± again approximates a correspondingly smaller constant term. Though the relative 
differences between these models are slightly more pronounced at larger scales, it is unlikely 
that current and emerging redshift surveys will attain sufficient precision to resolve them. 
Thus, while the geometric information contained in the present-day 3PCF can, unlike 5*3, be 
used to distinguish between different cosmological models and separate the contributions of 
constant linear- and nonlinear-bias terms, it may not reliably distinguish between constant 
and evolving bias. 

This ambiguity could, in principle, be removed by considering measurements of the 
spatial 3PCF as a function a redshift, as might be obtained from very deep surveys. Figure 
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Fig. 10. — Predicted Qv(v) f° r r = 8 h" 1 Mpc and u — 4 for the five bias scenarios 
shown. For each scenario, the ACDM and OCDM models are shown in thick and thin lines, 
respectively, with the two models being effectively identical in the non-evolving cases, but 
distinct in the evolving cases due to their different linear growth factors. 
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Fig. 11.— Predicted Q v (r,u,v) for the SCDM, 61 = 2.5, b 2 /b\ = 0.21 model (dot-dashed 
curves), OCDM, b u = 5, b 2 Jb 2 u = model (dashed curves), and OCDM/ACDM, 61 = 1.8, 
b 2 jb\ = model (dotted curves). The constant linear- and nonlinear- bias models (dotted 
and dot-dashed, respectively) yielded degenerate S3 predictions, but can be distinguished 
using the full geometric dependence of Qy (in part because of their different values of T). 
The evolving and non-evolving linear models (dashed and dotted, respectively), however, are 
not as well distinguished by measurements of the present-day 3PCF alone. 
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Fig. 12. — The variation of Qv with redshift for a fixed triangle configuration (r = 16 h^ 1 
Mpc, u — 1, v — 0.65). The dotted lines show the result for the OCDM/ACDM model with 
constant linear bias, b\ = 1.8, while the dashed lines show the evolving linear bias, &i„ = 5, 
result for the OCDM. The evolving and constant OCDM models are virtually identical at 
z = 0, but can be distinguished by their different redshift-dependences. An evolving, ACDM 
model is shown by the thin, solid lines, for comparison. The right panel shows the results 
obtained using a formation redshift of Zf = 10, rather than Zf = 5. 
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12 shows the evolution of Qy with redshift, for a particular configuration (r = 16 h" 1 Mpc, 
u = 1, v = 0.65), for the OCDM evolving linear bias model and OCDM/ACDM constant 
linear bias model from Figures |5] and [II], with a ACDM evolving linear bias model added 
for comparison. These evolving and non-evolving cases yielded similar results for Qy as 
well as 5*3, and for this particular configuration were practically identical, as can be seen in 



the appropriate panel of Figure |TI| or from the z = result in Figure 0. We see, however, 
that the constant linear bias model can be distinguished from the evolving models on the 
basis of their evolution with redshift. Note that, unlike the constant-bias case where the 
negligible dependence on \i caused the OCDM and ACDM predictions to differ only by the 
thickness of the curves, the evolving OCDM and ACDM models here yield different results, 
as expected, with the differences becoming larger for smaller values of Qq. Comparing the 
left and right panels further illustrates the dependence on the redshift of galaxy formation, 
Zf, which, in principle, can also be constrained. 

In theory, a survey large and deep enough to resolve this degeneracy between constant 
and evolving bias might also enable a mo del- dependent determination of b\ and 62 solely 
from measurements of S3, as described in §2.1. Since distinctions between evolving and 



non-evolving bias, such as those drawn in Figure [12], could also be made by considering 
the redshift dependence of S3, we can infer that such a survey might address the various 
degeneracies between b\, 62, and their evolution, solely from measurements of the skewness, 
without the need to consider the configuration dependence of the full 3PCF. Unfortunately, 
it is unlikely that upcoming redshift surveys will be both deep and large enough to resolve 
the degeneracy of ( with respect to constant vs. evolving bias, and in any case would still 
be susceptible to redshift distortions which mix the density and velocity fields. However, 
BKJ show that the angular 3PCF, which can be well-sampled over a large range of redshifts 
(e.g., by radio surveys) and is free of such distortions, might be better able to discriminate 
between evolving and non-evolving bias scenarios. 

Based on these results, it is clear that the often-stated result that Q is a constant of 
order unity (Groth & Peebles 1977; Jing & Zhang 1989; Toth, Hollosi, & Szalay 1989; Gott, 
Gao, & Park 1991; Jing, Mo, & Borner 1991; Baumgart & Fry 1991; Bouchet et al. 1993) is 
only a crude approximation resulting from a coarse averaging over allowed configurations 
and systematic errors arising from the normalization convention. Valuable information is 
lost if the geometric dependence of Q, or more appropriately Qy, is ignored. We also stress 
that Figures |9|-[l2|, are only illustrations; the statistical leverage gained by considering the 



full configuration-space variation of the 3PCF as a function of redshift should in principle 
discriminate between the signatures of gravitational evolution versus bias, independently 
measure the effective present-day values of b\ and 62, the shape of the power spectrum, 
and possibly measure the evolution in these quantities and constrain the epoch of galaxy 
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formation. 



4. CONCLUSION 

We have derived leading-order results for the normalized spatial skewness and 3PCF, 
assuming Gaussian ICs, for an arbitrary tracer-mass distribution in flat and open universes, 
taking into account such features as the scale dependence and linear evolution of the CDM 
power spectrum and the presence of a possibly evolving (decaying), linear or nonlinear bias, 
as defined for unsmoothed density fields. The predicted normalized amplitudes, S3 and 
Q, for an unbiased tracer mass are in agreement with previous work (Bouchet et al. 1992; 
Juszkiewicz, Bouchet, & Colombi 1993; FG93; Bernardeau 1994a; Fry 1994; JB97). We 
extend the FG93 result for S3 in the case of a nonlinear bias, defined for unsmoothed 
fields, to include a scale-dependent, leading-order correction which becomes appreciable 
for positive effective spectral indices, corresponding to scales R > 100 h' 1 Mpc for CDM 
models. This correction term implies that the value of the nonlinear bias parameter, as 
defined for smoothed density fields, could generally depend on the adopted smoothing scale. 
In the unlikely event that S3 could be measured over scales R > 100 h~ l Mpc, the presence 
or absence of this behavior could shed light on the smoothed vs. unsmoothed biasing 
picture, and perhaps allow more accurate determinations of the linear- and nonlinear-bias 
parameters on the basis of skewness measurements alone. 

We show that the conventional definition of Q gives rise to rapid, non-monotonic 
variation and divergences which do not arise from the behavior of the 3PCF. We speculate 
that the large discrepancy between this behavior and that predicted by some iV-body 
simulations is due to a practical bias, arising from the poor normalization convention, in the 
reconstruction of Q from data, rather than from the absence of nonlinear corrections in the 
QL PT prediction. We propose that this question may be addressed by considering iV-body 
results at earlier cosmic epochs. We consider instead Q V) which is everywhere well-behaved. 

We find the bias model to be a crucial factor influencing the clustering predictions, 
in agreement with Matarrese et al. (1997). In particular, the scale dependence of S3 and 
the configuration dependence of Qy bear characteristic imprints of bias, such as a relative 
flattening and decrease with increasing b\ and a relative increase with increasing b 2 - Unlike 
S 3 , the predictions for Q v are seen to depend significantly on the adopted CDM model, 
through the dependence on T, and for all models investigated, the dependence of the 3PCF 
on triangle geometry implies stronger clustering of smaller, elongated structures. The 
statistic S3, while preserving information about the overall scale dependence, is seen to 
provide only one test on the combination of the cosmological parameters, b\ (t), and b 2 (t), 
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yielding degeneracies between these quantities, particularly for smoothing scales R < 100 
h^ 1 Mpc. These are partially alleviated by considering the dependence of the full 3PCF 
on the triangle geometry. In particular, the variation of Qy can be used to distinguish 
between the effects of gravitational evolution and bias, place constraints on the value of T, 
and measure the effective constant bias terms b\ and hi- Evolving-bias models, however, 
are found to yield similar predictions for both Qy and S 3 as models with an appropriately 
smaller, constant bias, and thus cannot be reliably distinguished using spatial statistics at 
z = 0. Measurements of the 3PCF as a function of redshift could, in principle, directly 
measure bias evolution, Q , and the redshift of galaxy formation, but these might be 
more readily measurable using angular statistics (BKJ). A comparison of these predictions 
for S3 and Qy with several data sets characterizing different tracer populations would 
allow multiple, complementary constraints on the above parameters. The often-quoted 
"empirical" result, that the normalized 3PCF is a constant of order unity, is in part a result 
of coarse-averaging over all triangular configurations, which destroys much of this valuable 
information. 

In practice, other factors arise when comparing these predictions with observations. 
For example, redshift distortions, which mix information about the density and velocity 
distributions, must be accounted for when measuring spatial clustering (Fry & Gaztahaga 
1994; Jing & Borner 1998; Verde et al. 1998; Heavens, Matarrese, & Verde 1998; 
Scoccimarro, Couchman, & Frieman 1998). Finite- volume and boundary effects, as well as 
estimation biases, tend to reduce the observed clustering amplitude (Baugh, Gaztahaga, & 
Efstathiou 1995; Szapudi & Colombi 1996; Colombi, Szapudi, & Szalay 1998; Gaztahaga 
& Bernardeau 1998; Hui & Gaztahaga 1998), and Poisson noise (Gaztahaga 1994) and 
sampling variance must also be taken into account. In addition, the calculations above can 
be extended in various ways. Our assumed model for the bias is deterministic, but other 
stochastic models have been proposed which may yield different results, particularly on 
smaller scales (FG93; Catelan et al. 1998; Catelan, Matarrese, & Porciani 1998; Taruya, 
Koyama & Soda 1998). One might also investigate the effects of including higher-order 
nonlinear terms in the PT expansion, such as one-loop results (Jain & Bertschinger 
1994; Scoccimarro & Frieman 1996a,b; Scoccimarro et al. 1998); these are expected to be 
negligible on QL scales, but become increasingly important when comparing predictions 
with observations on scales where £ > 1. For example, iV-body simulations indicate that 
Q becomes relatively shape-independent in the nonlinear regime, suggesting that a linear 
bias is not the only mechanism capable of producing a flattening (Scoccimarro et al. 1998), 
and Gaztahaga & Bernardeau (1998) also find that nonlinear effects tend to increase S3 
while erasing the shape dependence of Q. Calculations analogous to those presented here 
can also be performed for higher-order moments and n-point CFs (e.g., kurtosis, four-point 
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correlation function, etc.) which, in principle, should provide additional constraints. Finally, 
the entire discussion can be extended to include the predictions of structure-formation 
models with non-Gaussian ICs, such as topological-defect or isocurvature models (e.g., Jaffe 
1994). 
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